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ABSTRACT 

The gas cloud G2 is currently being tidally disrupted by the Galactic Centre 
super-massive black hole, Sgr A*. The region around the black hole is populated by 
~ 30 Wolf-Rayet stars, which produce strong outflows. We explore the possibility that 
gas clumps, such as G2, originate from the collision of stellar winds via the non-linear 
thin shell instability. Following an analytical approach, we study the thermal evolution 
of slabs formed in the symmetric collision of winds, evaluating whether instabilities 
occur, and estimating possible clump masses. We find that the collision of relatively 
slow ( 750 km s -1 ) and strong (~ 10“ 5 M 0 yr -1 ) stellar winds from stars at short 

separations (< 10 mpe) is a process that indeed could produce clumps of G2’s mass 
and above. Such short separation encounters of single stars along their known orbits 
are not common in the Galactic Centre, making this process a possible but unlikely 
origin for G2. We also discuss clump formation in close binaries such as IRS 16SW 
and in asymmetric encounters as promising alternatives that deserve further numerical 
study. 
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1 INTRODUCTION 

Gillessen et al (2012) detected a moving diffuse object, 
the so-called G2 cloud, on its way towards Sgr A*, the 
radio source identified as the central massive black hole 
of our Galaxy (see Genzel et al. 2010, for a review). 
The tidal disruption of this cloud is being monitored by 
different groups (e.g., Eckart et al. 2013; Gillessen et al. 
2013; Phifer et al. 2013) and provides a unique opportunity 
to test accretion physics, due to both its proximity 
and the timescale on which it happens. Gillessen et al 
(2012) estimated the mass of G2 to be ~ 3M® from 
its line emission. However, the nature of G2 has not 
been clarified yet. In particular, there is an ongoing 
debate on whether the diffuse cloud contains a compact 
mass (likely a star). Witzel et al. (2014) presented the 
detection of a compact source at 3.8 fim (thermal dust 
emission) that would correspond to G2 during pericentre 
passage. Its survival as a compact source to the close 
passage (~ 2000 Schwarzchild radii) suggests the existence 
of a central star keeping it bound, which they argue 
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is a binary merger product (Prodan et al. 2015). Other 
possible explanations for a central mass in G2 include an 
evaporating proto-planetary disc or the wind of a T-Tauri 
star (Murray-Clay & Loeb 2012; Scoville & Burkert 2013; 
Ballone et al. 2013; Do Colle et al. 2014). 

On the other hand, Pfuhl et al. (2015) argued in favour 
of a purely gaseous cloud nature for G2 using their 
Brackett-y observations. They interpreted this source as 
a bright knot of a larger gas streamer that includes a 
G2-type object called G1 in a similar orbit, but preceding 
it by 13 yr. G1 and G2 could be explained as the result 
of the partial tidal disruption of a star (Guillochon et al. 
2014) or as one of many gas clumps created by the 
collision of stellar winds from the young stars in the 
Galactic centre (Burkert et al. 2012; Schartmann et al. 
2012). Such dense, cold clumps are copiously produced 
in the Smoothed-Particle Hydrodynamics (hereafter SPH) 
simulations of the Galactic centre gas dynamics performed 
by Cuadra et al. (2005, 2006, 2008, 2015) (see also 
Liitzgendorf et al. 2015 § 3.3), and could survive pericentre 
passage if magnetised (McCourt et al. 2015). Moreover, G2’s 
orbit lies on the plane of the ‘clockwise disc’, defined 
by the orbits of many young stars (Paumard et al. 2006; 
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Yelda et al. 2014), and its apocentre coincides with the inner 
rim of that disc. Nevertheless, the SPH technique has a 
tendency for artificially clumping gas (Hobbs et al. 2013), 
which raises doubts on how physical the clump formation 
is in such models. In this context, this work aims to test 
independently the clump formation as result of colliding 
stellar winds in the central parsec of the Milky Way. 

Massive stars have significant phases of mass loss 
through their lives, in which their outflows can be 
accelerated up to supersonic speeds due to radiation 
pressure. When two of these stars are at short separations, 
for example in a binary system, their winds collide and 
generate a hot slab of shocked gas. Depending on the ability 
of this gas to cool down, we can identify different regimes. If 
the shocked gas from both stellar winds cool down rapidly, 
a cold thin shell will be produced centred at the contact 
discontinuity (hereafter CD) of the wind-wind interaction 
zone. This slab will be subject to strong instabilities such 
as the Non-linear Thin Shell Instability (hereafter NTSI, 
Vishniac 1994). In the case that only one of the winds 
is highly radiative, a cold thin shell will also be formed; 
but if instabilities are excited they will be damped by the 
thermal pressure of the hot shocked gas from the other 
wind (Vishniac 1983). When none of the winds are radiative, 
the Kelvin-Helmholtz instability (hereafter KHI) is the only 
instability that can be excited (Stevens et al. 1992), as 
it only requires a velocity difference between the winds. 
However, KHI can also act on top of thin-shell instabilities 
for radiatively efficient winds. 

State-of-the-art numerical modelling (Pittard 2009; 
van Marie et al. 2010; Lamberts et al. 2011; Kee et al. 
2014) highlights the high computational cost of realistic 
simulations of unstable colliding wind systems. High spatial 
and time resolution, plus many physical ingredients, such as 
gravity, driving of the winds, radiative cooling and orbital 
motion, are crucial to build realistic models and track the 
growth of instabilities. In this context, we take an alternative 
analytical approach and study the thermal evolution of the 
hot slab created by the colliding winds. This allows us to 
predict under which conditions, out of a wide parameter 
space, the NTSI grows; and to estimate the possible resulting 
clump masses. The paper is divided as follows: Section 2 
describes the main cooling diagnostic we use throughout 
this work and Section 3 presents our model and the results 
of our parameter space study. In Section 4 we discuss how 
likely the formation of the G2 cloud is through colliding 
winds based on our results. Then in Section 5 we discuss 
the limitations of our results and argue for colliding winds 
binaries as clump sources in the Galactic Centre. Finally, in 
Section 6 we present our conclusions and outlook. 


2 RADIATIVE COOLING AND THIN SHELLS 

The role of radiative cooling is what determines the 
thickness, density and temperature of the shocked gas 
layer. If cooling is efficient, the wind collision will produce 
a dense, thin layer of cold gas, which can easily be 
subject to the so-called thin shell instabilities. In their 
numerical models, Stevens et al. (1992) identified two such 
instabilities: the NTSI, later explained analytically by 
Vishniac (1994), and another damped instability consistent 


with the earlier description by Vishniac (1983). Moreover, 
Dgani et al. (1993) described another thin shell instability: 
the Transverse Acceleration Instability (TAI) that mainly 
takes place off the two-star axis. Lamberts et al. (2011), 
modelling unstable colliding winds systems numerically, 
concluded that the NTSI is the instability that dominates 
the cold slab evolution due to its large scale perturbations. 

In order to predict whether the wind collision will 
produce a thin slab we need to check how fast the shocked 
gas cools down. Following the description presented in 
Stevens et al. (1992), the radiative efficiency of shocked 
stellar winds can be described by the cooling parameter X- 
which is the ratio of the cooling timescale to the dynamical 
timescale, 


X = 


f cool 
tdyn 


( 1 ) 


where t coo i and ta yn are the cooling and the dynamical 
timescales, respectively. They are defined as, 


3ksT d* 

tco ° l = 2nA(T) ; Uyn = ^ 


( 2 ) 


where T and n are the temperature and number density of 
the shocked gas, respectively, A(T) is the cooling function 
which we will specify in the following section, fcn is the 
Boltzmann constant, d* is the distance between the star and 
the CD and c s is the shocked gas sound speed. Therefore, 
for x lower than unity, the gas can cool down faster through 
radiative cooling than through adiabatic expansion. On the 
contrary, \ larger than one means that radiative cooling 
is not efficient. Following Stevens et al. (1992), but using a 
cooling function that depends on metal abundance (specified 
in Section 3.2) and temperature rather than being constant 
(see Equation 10) we can calculate the cooling parameter 
for each stellar wind given its properties, 

1 V s 5 ' 4 d *12 


X ~ 2 " 


M- 7 


( 3 ) 


where Vs is the terminal wind speed V in units of 
1000 km s -1 , d* 12 is the distance between the star and the 
CD d* in units of 10 12 cm and M_ 7 is the mass loss rate 
M in units of 10 -7 M@ yr -1 . In Figure 1, we show x as 
a function of V and the stellar separation of identical stars 
d = 2d* for M = 10 -5 Mq yr -1 fixed, which is the typical 
value for the Wolf-Rayet stars in the Galactic centre (see 
§ 4). We have highlighted the boundary x = 1 to separate 
the two cooling regimes. We will concentrate on the x < 1 
region, which is where the NTSI excitation is possible. 


3 NTSI AND CLUMP FORMATION 

As discussed in the previous section, if the slab gas can 
cool down rapidly, different thin shell instabilities can be 
excited. The dominant instability is the NTSI, which is the 
result of the misalignment of the thermal pressure within 
the cold slab (which always acts perpendicular to the slab) 
and the ram pressure of the wind (always acting parallel to 
the wind direction). This generates a convergent flux onto 
perturbation knots where gas accumulates, and a subsequent 
mixing of both phases of the material. In this section we use 
the description of Vishniac (1994) to calculate the allowed 
wavelength range of the NTSI which can grow. 
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Figure 1. Cooling parameter as colour map on logarithmic scale 
which was obtained using Equation 3 and fixing the mass loss 
rate to 10 -5 Mq yr" 1 . It is shown as a function of the stellar 
separation d (for identical stars) on the x-axis and the stellar 
wind velocity V on the y-axis. The solid line shows X = 1 which 
separates the two regimes: adiabatic (above, x > 1) and roughly 
isothermal (below, \ < 1) winds. Dashed and dotted lines show 
the x = 1 boundaries for mass loss rates of 10“ 4 and 10“ 6 
Mq yr -4 , respectively. Green stars labelled A-H indicate the 
parameters of models we study in more detail in § 3. Notice that 
1 mpc ~ 206 AU. 


3.1 Range of unstable wavelengths 


slab. All these quantities evolve and we study them as a 
function of their age, i.e., the time after the wind collision. 
As we are studying a symmetric system (identical stars with 
identical winds), we model one side of the system. Therefore, 
L(t) will be the distance from the CD to the discontinuity 
of the shocked and the free wind region. As the slab is 
supported by thermal pressure, in this point the thermal 
and the wind ram pressure have to be balanced, which 
is a reasonable assumption considering that the thermal 
energy is negligible in the free-wind region that is highly 
supersonic (M ^ 20). From this model, we can write the 
hydrodynamical equations as follows, 


4.5 

7 , S s lab 
dt 

— Pwind^4 

( 5 ) 

6.0 

Psiab 

= Pwind ^ ? 

(6) 


4^-4 (SslabT s l ab ) 

Ifxm h dt 

= -^shock + *Scool} 

( 7 ) 


where P s iab is the pressure 

in the slab, p is the 

mean 


molecular weight, wh is the proton mass, ffshock is the 
mechanical heating term (i.e. kinetic energy flux from the 
wind per unit surface), and S coo \ is the energy dissipation 
term through radiative cooling. These source terms are given 
by 


7 7shock — 


'Scool — 


, PwindE , 

E s lab Psiab 

Cml 


A(T’ a lab), 


( 8 ) 

( 9 ) 


where the cooling function A (T) is an analytical 
approximation for optically thin radiative cooling made by 
Cuadra et al. (2005) following Sutherland & Dopita (1993), 


The growth of the NTSI can occur only for wavelengths 
shorter than those coordinated by sound waves: A < c s t. 
On the other hand, the shortest unstable wavelength is 
given by the slab thickness l(t), i.e., A 2/(f). Moreover, 
the NTSI can grow only if the slab perturbation has an 
initial amplitude at least comparable to the thickness of the 
layer. Notice that the criteria depend on the time t elapsed 
since the formation of the slab. The thickness of the slab 
generated by isothermal flows is given by l(t) = 2V n t, where 
V n = V/{M 2 — 1) is the velocity at which the slab thickens, 
Ad = Va/ca is the Mach number, and 14 = V + V n is the 
shock speed. Notice that for the limit of high Mach number 
Vs ~ V, and that the more supersonic the shock, the slower 
the slab increases its width. 

3.2 Analytic model 

In our model we consider two identical stars that are fixed in 
space. They have a given mass loss rate and wind terminal 
velocity, and are separated by a given distance. Respectively, 
these three quantities M, V and d, are the input parameters 
of the model. We follow the time evolution of the gas 
within the slab formed when the winds collide by studying 
its integrated density and thermal evolution, including a 
radiative cooling term. We define the relevant densities as, 

E s lab — PslabT, (4) 

where p a i a b is the volumetric density inside the slab, E s i a b 
is the surface density of the slab and L is the width of the 


A (T) = 6.0 x 10- 23 (|) 0 7 erg cm 3 s" 1 , (10) 

where Z is the metal abundance relative to solar. We set 
Z = 3 and discuss the impact of this choice in Section 5.1.2. 
Combining the system of Equations 5, 6, 7, 8, 9; and 
assuming an ideal gas P s iab = Psiabfc B Tsiab/(pm H ), we 
can derive a single differential equation that describes the 
thermal evolution of the slab: 

3fc B d ^ ^ PwindC 2 A(T’ s lab) 3 fc B Tslab 

2/rrriH dt 2 t fc B pmH T s i a b 2pm,H t 

( 11 ) 

The first term on the right hand side represents the 
mechanical heating of the shock; the second, the radiative 
cooling; and the third, the work done by the slab. 
Furthermore, the slab density can be easily computed 
making use of Equation 6 and the ideal gas equation of state. 
Thus, the slab width time evolution can be calculated using 
Equations 4 and 5 to obtain 

4 [Psiab (t)L(t)} = p wind c => L(t) = -£2±rVt. (12) 

dt p s l a b ( t) 

From this expression we note that if the slab does not 
cool down (i.e. the slab density remains roughly constant), 
the slab will increase its width linearly with time; on 
the contrary, if the slab gets denser very rapidly it may 
overcome the linearity dependence with time and collapse 
as a thin shell which grows very slowly due to the fact that 
Pwind/pslab <S 1- With this analysis we can calculate the 
NTSI growth criteria as a function of the age of the slab (or 


MNRAS 000, 1-11 (2015) 

















4 D. Calderon et al. 


Table 1. Parameters of each model presented in this work. 
Column 1: Model name. Column 2: wind terminal velocity. 
Columns 3 and 4: stellar separation in mpc and AU, respectively. 
Column 5: cooling parameter y estimated from Equation 3 and 
M fixed to 10 -5 Mg yr -1 . 


Model 

V 

(km s -1 ) 

(mpc) 

d 

(AU) 

X 

A 

250 

10.0 

2.06 x 10 3 

0.0433 

B 

500 

1.0 

2.06 x 10 2 

0.1827 

C 

250 

1.0 

2.06 x 10 2 

0.0043 

D 

750 

0.1 

2.06 x 10 1 

0.1632 

E 

500 

0.1 

2.06 x 10 1 

0.0183 

F 

250 

0.1 

2.06 x 10 1 

0.0004 

G 

1500 

1.0 

2.06 x 10 2 

69.89 

H 

2000 

0.1 

2.06 x 10 1 

32.57 


time after the wind collision). Although the model is simple, 
it can provide us with information about the instabilities 
that can take place in order to estimate sizes and masses 
of the possible clumps. It is important to mention that 
throughout this analysis we are not attempting to follow 
the evolution of the instability, but to study under which 
conditions it can be triggered. 

3.3 Parameter study 

Using the previously described procedure, we explored the 
parameter space shown in Figure 1. We first present in detail 
a few relevant and representative examples (see Table 1). 
In all these models we kept M = 10~ 5 Mq yr -1 fixed 
and stopped the calculation after the gas has cooled down 
significantly, reaching a temperature of 10 4 K. This is a 
reasonable temperature floor because of the presence of 
many hot, young stars emitting ionising radiation that 
prevents the temperature to drop to lower values 1 . 

The time evolution of the unstable wavelength criteria 
for the NTSI are shown for all models in Figures 2 and 3. 
Figure 2 shows the models with y < 1. Here, for most of 
the evolution of the slab, the upper limit (solid blue line) is 
below the lower limit (dashed green line), meaning that the 
instability cannot grow as there is no unstable wavelength 
range. Only at the very end of each calculation the upper 
limit is above the lower limit and therefore the instability 
develops. The sudden drop in the lower limit is due to 
the fast slab gas compression once cooling becomes much 
more efficient as the temperature decreases. In all these 
cases we registered the unstable length range and the final 
density values in Table 2. It is important to remark that 
for all these cases we remain in the thin shell regime, i.e., 
L/d <C 1. Therefore we do not expect radiative cooling to 
be faster than the thermal response of the slab, making our 
assumption of pressure equilibrium at the shock valid. 

Assuming spherical symmetry, we used the unstable 
wavelength range to compute a range of masses for clumps 
(see Table 2). However, these values must be interpreted 

1 We also tested stopping the cooling at 10 s K, finding 
that density values are systematically smaller by one order of 
magnitude. 


only as upper limits, because we have assumed the clump 
radii to be about the size of the unstable wavelengths, 
which is probably an extreme case. From the results we see 
that larger stellar separations will result in larger clumps. 
However, these clumps are less dense due to the winds being 
significantly diluted before they collide. On the contrary, 
smaller separations produce smaller and denser clumps. The 
combination of these two factors sets the clump masses in a 
non-trivial manner. Our results show that Model A, Model B 
and Model D could generate clumps with masses > 0.1 Mq 
and the most massive ones would be generated by Model B, 
reaching G2-like masses (~ 5 Mq). 2 

On the other hand, Figure 3 shows the models with y > 
1 , in which we see that the width of the slab (dashed green 
line) becomes larger than the stellar separation (dotted red 
line) before there is an unstable wavelength range. Thus, we 
deem these models not physical, as our treatment cannot 
describe these systems properly, and we expect no clump 
formation through NTSI for this parameter range. 

Although these results show that clump formation might 
occur for the y < 1 models, we need to check on which 
timescales this process takes place. Vishniac (1994) showed 
that the NTSI growth timescale r is given by 


where C, is the initial amplitude of the perturbation. If we 
assume £ ~ A, the growth timescale is simply given by the 
sound crossing timescale, i.e., r ~ A/c s . Estimations of r 
under this assumption for each of our models are presented 
in Table 2. Moreover, we present the clump masses as a 
function of the growth timescale for each model in Figure 4, 
where we see that models with shorter stellar separations 
(d = 0.1 mpc in green lines) tend to be the ones that can 
create clumps the quickest. This is due to these winds being 
less diluted when they collide, making the cooling more 
efficient. We can compensate this effect by increasing the 
wind speed, as a hotter slab would take longer to radiate 
most of its energy away. For example, comparing cases C 
and E (solid blue and dashed green, respectively) we see that 
the combination of different stellar separations and different 
wind speeds produces clump formation on roughly the same 
timescale. 

An interesting point is that the models that can 
create the most massive clumps have significantly different 
parameters. This can be explained by the fact that clump 
masses are proportional to A 3 , thus the most massive clump 
possibly generated would have a mass M oc A^ ax . Moreover, 
the upper limit of the unstable wavelength range of the NTSI 
is given by the sound crossing distance, i.e., f c s dt integrated 
over the age of the slab when it reaches 10 4 K. This tells 
us that how massive clumps can be, depends on how long 
it takes for the slab to cool down and become unstable. 
Then, either high velocity winds and/or larger separations 
would produce more massive clumps. In order to illustrate 
this explicitly we have explored our parameter space more 
extensively, modelling 45 systems in total. We show this in 
Figure 5 as a function of the input parameters of our model. 

2 G2 could also have formed by mergers of smaller clumps, but 
we cannot address this option with our current approach. 
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Model A : V — 250 kin s 1 and d — 10.0 mpc => \ — 0.0433 



Model C : V — 250 kin s 1 and eZ — 1.0 mpc => y — 0.0043 



Model E : V = 500 km s 1 and d = 0.1 mpc =► x = 0.0183 
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Model B : V = 500 km s 1 and d = 1.0 mpc ;■ y — 0.1827 



Model D : V = 750 km s 1 and d = 0.1 mpc y — 0.1632 



Model F : V = 250 km s 1 and d = 0.1 mpc => x = 0.0004 
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Figure 2. Evolution of the NTSI wavelength criteria for the radiative models (y < 1). On the upper panels we show the time evolution 
of the upper and lower limits of unstable wavelengths for the NTSI obtained from our semi-analytic prescription. Instabilities can grow 
only at the end of each time evolution, when a range of allowed wavelengths exists. Yellow stars represent the final state of the system 
when the gas has reached 10 4 K. On the lower panels, the thick black line represents the difference between the upper and lower limits, 
A max — A min , and the dotted red line is fixed at zero for reference. Note the different scales in the x- and y-axis between models. 


The left and right panels present the lower and upper limits 
of the clump mass range, respectively. Here it is easy to see 
that for larger separations and/or higher wind velocities the 
minimum and maximum clump masses are larger. 

We repeated the previous analysis using different mass 


loss rate values, 1CP 6 and 10~ 4 Mq yr -1 , finding that for 
lower M values, clumps can be more massive for the same 
combination of ( d, V), as a less dense slab makes the cooling 
less efficient, but we would need shorter distances and/or 
slower winds in order to get radiative winds in the first 
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Model G : V — 1500 kin s 1 and d — 1.0 mpc => x~ 68.89 



Slab Age t (yr) 


Model H : V — 2000 kin s 1 and d — 0.1 mpc => \ — 32.57 



Figure 3. Same as Figure 2 but for the adiabatic models (x > 1) G and H. At some point the models become unphysical, as the slab 
width is comparable to the stellar separation (red dotted line). We do not expect NTSI excitation for these models. 


Table 2. Results of our estimates for each model. Column 1: model name. Columns 2 and 3: unstable wavelength range obtained from 
our analysis in mpc and AU, respectively. Column 4: slab density obtained estimated from our calculations. Column 5: clump mass range 
computed assuming spherical clumps and uniform density with radius equal to the instability wavelengths. Column 6: growth timescale 
estimated assuming an initial amplitude equal to the instability wavelength and the slab sound speed at the temperature floor (10 4 K). 


Model 

Unstable A 

Pslab 

Clump Mass Range 

Growth Timescale 


(mpc) 

(AU) 

(g cm -3 ) 

(M e ) 

(yr) 

A 

(4 - 190) X 10“ 3 

(8 - 390) x 10" 1 

2.4 X 10 -18 

3.6 x 10“ 6 - 3.3 x 10 _1 

(3 - 120) x 10 _1 

B 

(5 - 800) x 10~ 4 

(1 - 160) x 10" 1 

4.7 x 10~ 16 

8.8 x 10 -7 -4.9 x 10° 

(3 - 50) x lO^ 1 

C 

(4 - 190) x 10 -5 

(8 - 390) x 10" 3 

2.4 x 10~ 16 

3.6 x 10 _1 ° - 3.3 x 10~ 5 

(3 - 120) x 10“ 3 

D 

(2 - 710) x 10 -5 

(4 - 1460) x 10 -3 

7.1 x 10" 14 

8.5 x 10“ 9 - 5.3 x 10 _1 

(1 - 470) x 10“ 3 

E 

(5 - 800) x 10~ 6 

(1 - 160) x 10“ 3 

4.7 x 10" 14 

8.8 x 10~ n - 4.7 x 10 -4 

(3 - 50) x 10" 4 

F 

(4 - 190) x 10 -7 

(8 - 390) x 10" 5 

2.4 x 10~ 14 

3.6 x 10~ 14 - 3.3 x 10~ 9 

(3 - 120) x 10“ 5 

G 

- 

- 

2.6 x 10" 19 

- 

- 

H 

- 

- 

2.6 x 10 -21 

- 

- 


place, as the x = 1 line moves down in the parameter space 
(see Figure 1). On the contrary, increasing M we obtain 
less massive clumps for the same (d, V), as a denser slab 
makes the cooling more efficient, but this would allow larger 
separations and/or faster winds to generate radiative winds, 
as the x = 1 line moves up in the parameter space. 

To study which clumps will actually form, we take into 
account Equation 13 that shows that shorter wavelengths 
grow faster. Therefore for any model we would expect first 
the formation of the lightest possible clumps. On longer 
timescales, larger wavelengths would also act, accumulating 
the small clumps and possibly merging them to create 
more massive clumps. To actually predict a clump mass 
distribution, we require numerical simulations. We defer 
that study to a forthcoming work. 

In general, we find that for parameters closer to x = 1 
(black dashed line in Figure 5) we would expect more 
massive clumps to be formed. However we do not know 
exactly at what point our approach becomes unphysical as 
we go closer to the adiabatic regime, as seen in models G 
and H (see Figure 3). Thus, it is more sensible to explore the 
parameter space close to \ = 1 with numerical simulations. 
This is why we explored values up to x = 0.5 only (black 
solid line). Despite this, we see that clump masses span a 


very wide range of masses and that the creation of clumps 
as massive as G2 is possible. 


4 COLLIDING STELLAR WINDS IN THE 
GALACTIC CENTRE 

We have studied the parameter space for close pairs of 
mass-losing stars, checking whether we expect them to 
form clumps through the NTSI, plus estimating the allowed 
range of sizes and masses for those clumps. Now, we apply 
this model to the stars in the central parsec of the Milky 
Way, as the NTSI could be excited for the colliding winds 
of massive stars, and explain the origin of the G2 cloud. 
Our model assumes that the stars are stationary and that 
the pairs are identical. The first assumption is justified as 
the time-scales involved in the instability (,$10yr, e.g., 
Figure 4) are typically much shorter than the duration of 
close stellar encounters of mass-losing stars in the Galactic 
Centre (see Table 4). The second issue is partially addressed 
in Section 4.3. Our sample of mass-losing stars in the inner 
parsec of the Galaxy is the same used in the numerical 
models by Cuadra et al. (2008, 2015) and it is listed in Table 
3. For most of the stars, the wind properties were taken 
from Martins et al. (2007), who fitted stellar atmosphere 
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Figure 4. Clump masses (assuming spherical symmetry) 
generated through the NTSI, as a function of the instability 
growth timescale (assuming an initial amplitude comparable 
to the instability wavelength). Different stellar separations are 
shown with different colours: green, blue and black stand for 0.1, 1 
and 10 mpc (20.6, 206 and 2060 AU), respectively, while different 
wind velocities are shown with different line styles: solid, dashed 
and dotted represent 250, 500 and 750 km s — . 


models to spectra obtained with SINFONI at the VLT. For 
the rest, the wind properties were assigned by Cuadra et al. 
(2008) based on their similarity to other stars whose spectra 
were properly modelled by Martins et al. (2007). The typical 
uncertainty in the parameters derived by Martins et al. 
(2007) is ~ 100 km s -1 for the velocities and ~ 0.2 dex 
for the mass-loss rates, small enough not to influence our 
conclusions. 


4.1 Radiative cooling diagnostic 

We first estimate the critical separation between stars in 
order for their winds to be radiatively efficient. We use the 
previously defined cooling parameter y, and compute the 
radiative-to-adiabatic transition separation equating y = 1, 

.cool n , i 7 ,\ 

d» =2x10 cm—. (14) 

V 8' 

Pairs of identical stars with separations d > 2d*° o1 will 
produce adiabatic shocks, while those with d < 2d* ot>1 will 
result in radiatively-cooled shocks. The value of 2d*° o1 is 
included in Table 3, and is also plotted as a function of the 
wind momentum fluxes in Figure 6. Note that, for any pair 
of stars, d changes with time as they orbit in the Galactic 
Centre, therefore colliding wind stars can go through both 
regimes. Star 33E has the largest value of d»° o1 . Encounters 
involving this star with similar ones at 10 mpc-scales should 
result in the formation of cold slabs. For the rest of the 
sample we see that only for separations below 0.1 — 1 mpc 
(20-206 AU) their winds would be radiatively efficient. Such 
short separations are in the range of close binary systems 
(Section 5.2). 


Table 3. Galactic Centre mass-losing star sample taken from 
Martins et al. (2007) and Cuadra et al. (2008). Column 1: star 
ID. Column 2: star name. Both from Paumard et al. (2006). 
Column 3: stellar wind terminal velocity. Column 4: stellar mass 
loss rate. Columns 5 and 6: twice the radiative-to-adiabatic wind 
transition distance computed from Equation 14 in mpc and AU, 
respectively. 


ID 

Name 

V 

M 

2 d c ° o1 



(km s _1 ) 

(M 0 yr- 1 ) 

(mpc) 

(AU) 

19 

16NW 

600 

1.12 X 10“ 5 

2.29 

472 

20 

16C 

650 

2.24 X 10 -5 

2.97 

612 

23 

16SW 

600 

1.12 X 10 -5 

2.29 

472 

31 

29N 

1000 

1.13 X 10“ 5 

0.15 

31 

32 

16SE1 

1000 

1.13 X 10“ 5 

0.15 

31 

35 

29NE1 

1000 

1.13 x 10 -5 

0.15 

31 

39 

16NE 

650 

2.24 x 10 -5 

2.97 

612 

40 

16SE2 

2500 

7.08 X 10“ 5 

0.01 

2 

41 

33E 

450 

1.58 X 10“ 5 

15.28 

3148 

48 

13E4 

2200 

5.01 x 10 -5 

0.01 

2 

51 

13E2 

750 

4.47 x 10 -5 

2.74 

564 

56 

34W 

650 

1.32 x 10“ 5 

1.75 

361 

59 

7SE 

1000 

1.26 X 10“ 5 

0.16 

33 

60 

- 

750 

5.01 x 10“ 6 

0.31 

64 

61 

34NW 

750 

5.01 x 10 -6 

0.31 

64 

65 

9W 

1100 

4.47 X 10 -5 

0.35 

72 

66 

7SW 

900 

2.00 X 10“ 5 

0.46 

95 

68 

7W 

1000 

1.00 x 10 -5 

0.13 

27 

70 

7E2 

900 

1.58 x 10 -5 

0.36 

74 

71 

- 

1000 

1.13 X 10“ 5 

0.15 

31 

72 

- 

1000 

1.13 X 10“ 5 

0.15 

31 

74 

AFNW 

800 

3.16 x 10 -5 

1.37 

282 

76 

9SW 

1000 

1.13 x 10 -5 

0.15 

31 

78 

B1 

1000 

1.13 X 10“ 5 

0.15 

31 

79 

AF 

700 

1.78 X 10“ 5 

1.58 

325 

80 

9SE 

1000 

1.13 x 10 -5 

0.15 

31 

81 

AFNWNW 

1800 

1.12 x 10" 4 

0.06 

12 

82 

Blum 

1000 

1.13 X 10“ 5 

0.15 

31 

83 

15SW 

900 

1.58 X 10“ 5 

0.36 

74 

88 

15NE 

800 

2.00 x 10 -5 

0.87 

179 


4.2 Mass-losing stars encounters in the Galactic 
Centre 

Now we turn our attention to the stellar orbits in the 
Galactic centre, to study how often close encounters between 
single mass-losing stars are produced . To do so, we ran a 
simple, Newtonian gravity, test-particle simulation to follow 
the stars around Sgr A*. This approximation is correct, as 
the distances from the stars to Sgr A* are too large for 
relativistic effects to be important, and the time-scale for 
stellar scattering to be relevant is much longer than the 
period we are interested in (see Alexander 2005). We ran the 
models for 10 4 yr, which is already a much longer timescale 
than the expected lifetime of a cloud like G2 (Burkert et al. 
2012 ). 

We only registered encounters of stellar pairs at distances 
shorter than 10 mpc (2060 AU) based on our parameter 
space study, noting the minimum distance of the close 
passage d m i n and the duration of the encounter (i.e. the 
time stars are closer than 10 mpc). As initial conditions 
we used the 3D velocities and 2D sky positions observed 
by Paumard et al. (2006), meanwhile the unobservable 
2 -coordinate was chosen using different assumptions for the 
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M = lO- 5 Mguu yr- 1 



Stellar Separation d (mpc) Stellar Separation d (mpc) 


Figure 5. Parameter space ( d , V) with M = 10 —5 Mq yr —1 fixed, displaying clump masses (colours) formed through NTSI obtained with 
our model. The left and right panels show the minimum (lower limit) and maximum (upper limit) clump mass for a given combination 
of parameters, respectively. Green star symbols show the models we previously studied in detail. The solid and dashed black lines stand 
for x = 0.5 and x — 1> respectively. 


10 ' 


10 1 


w 


I* 10- 1 


10 “ 


10 “ 


★ 

33E 



vV 

★ 



V 


★ 

AFNWNM 


★ 

13E4 

★ 

16SE2 


10 4 


10 3 


10 2 < 


10 1 


10 ° 


0.00 0.05 0.10 0.15 0.20 0.25 

Momentum Flux (Msun yr -1 km s 1 ) 


Figure 6. Critical stellar separation in order for their winds to 
be efficiently radiative (i.e. their collision could produce a thin 
shell ) as a function of their wind momentum flux (M X V). 
The sample plotted corresponds to the one shown in Table 3. 
Stars with large momentum fluxes are labelled, as they can be 
important in asymmetric encounters “forcing” weaker winds to 
radiate their energy rapidly, resulting in radiative shocks. Star 
33E is also labelled as it has the largest value of d“° o1 . 


orbital distribution. We used the 3 different models from 
the work by Cuadra et al. (2008) to obtain the ^-coordinate 
of our young star sample: Min-ecc, IDisc and 2Discs. 
Min-ecc uses the ^-coordinate values that minimise the 
orbital eccentricities, IDisc assumes roughly half of the stars 
are in the well defined clockwise disc (Beloborodov et al. 
2006), and 2Discs assumes the existence of both clockwise 
and counter-clockwise discs. The results of this procedure 
are summarised in Table 4, where we see that over a 


10 4 yr period, there is at most one encounter with a ~ 1 
mpc separation 3 . From the typical encounter duration of 
~ 50 yr and our integration time, we can estimate a rough 
probability of producing G2 through encounters of single 
stars as 0.5%. 


4.3 Asymmetric mass losing stars encounters 


For simplicity, we have only analysed collisions of identical 
stellar winds. However, in reality we will typically have 
encounters between stars with different wind properties. To 
study such cases, detailed numerical simulations are needed, 
such as those performed by Pittard (2009); van Marie et al. 
(2010); Lamberts et al. (2011). These authors have shown 
that even a small velocity difference in the colliding winds 
can excite the KHI. This would mix two-phase material on 
top of other instabilities that can take place simultaneously. 
Although all these processes can be very complicated 
to track analytically, based on our work we can give 
a qualitative description of possible scenarios that can 
take place for asymmetric close encounters. Stellar wind 
encounters are characterised by their momentum flux ratio 
(Lebedev & Myasnikov 1990), 


M2V2 

M 1 V 1 ’ 


(15) 


where the subscript 2 stands for the weaker wind and 1 for 
the stronger one, so 77 ^ 1 by definition. Notice that 77 is 
independent of the stellar separation, provided this is large 
enough to allow the winds to reach their terminal velocities. 
For 77 ^ 1, the interaction zone of the shocked gas bends 


3 Notice that this “encounter” happens between stars 13E2 and 
13E4, which might be bound together by a dark mass (Fritz et al. 
2010), not included in our calculations. That dark mass would, 
however, likely increase the encounter duration. 
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Table 4. Galactic Centre mass-losing star encounters at < 10 mpc (2060 AU). Column 1: initial conditions model name. Column 2: IDs 
of the stars in the encounter. For each pair, the asterisk symbol (*) marks the star with the weaker wind (i.e., smaller momentum flux). 
When no asterisk is shown, both stars have the same momentum flux. Columns 3 and 4: minimum stellar separation in the encounter 
in mpc and AU, respectively. Column 5: time stars are closer than 10 mpc (2060 AU). Column 6: wind momentum flux ratio estimated 
from Equation 15. Columns 7 and 8: CD distance to the weaker wind star calculated from Equation 16 in mpc and AU, respectively. 
Column 9: weaker wind cooling parameter obtained from Equation 17. Column 10: whether the weaker wind is radiative or not at the 
minimum stellar separation of the encounter. 


Model 

Stars 

^min 

(mpc) (AU) 

Duration 

(yr) 

V 

R 2 

(mpc) 

(AU) 

Weaker wind 

Xasym 

Thin shell? 
(Xasym < 1) 

IDisc 


40 - 60* 

9.2 

1895 

10 

0.021 

1.16 

239 

7.581 

NO 


48 - 60* 

4.0 

824 

80 

0.034 

0.62 

128 

4.052 

NO 


48 - 51* 

1.0 

206 

1200 

0.304 

0.36 

74 

0.260 

YES 


51 - 60* 

6.0 

1236 

70 

0.112 

1.50 

309 

9.797 

NO 

2Discs 


19* - 48 

3.0 

618 

15 

0.061 

0.59 

122 

0.519 

YES 


19* - 51 

7.0 

1442 

20 

0.200 

2.16 

445 

1.889 

NO 


19* - 48 

2.0 

412 

220 

0.061 

0.40 

82 

0.347 

YES 


19* - 51 

4.0 

824 

20 

0.200 

1.24 

255 

1.080 

NO 


19* - 51 

9.5 

1957 

10 

0.200 

2.94 

606 

2.564 

NO 


20 - 32* 

8.0 

1648 

60 

0.776 

3.75 

773 

51.16 

NO 


31 - 72 

7.5 

1545 

20 

1.000 

3.75 

773 

51.20 

NO 


39* - 51 

2.0 

412 

15 

0.434 

0.79 

163 

0.534 

YES 


41* - 48 

9.0 

1854 

10 

0.065 

1.83 

377 

0.240 

YES 


48 - 51* 

5.0 

1030 

550 

0.304 

1.78 

367 

1.297 

NO 

Min-ecc 


23* - 41 

8.0 

1648 

50 

0.945 

3.94 

812 

3.442 

NO 


48 - 51* 

1.0 

206 

1200 

0.304 

0.36 

74 

0.260 

YES 


48 - 60* 

4.0 

824 

80 

0.034 

0.62 

128 

4.058 

NO 


48 - 61* 

6.0 

1236 

20 

0.034 

0.93 

192 

6.084 

NO 


51 - 60* 

6.0 

1236 

70 

0.112 

1.50 

309 

9.797 

NO 


51 - 61* 

7.5 

1545 

25 

0.112 

1.88 

387 

12.25 

NO 


towards the weaker star. The CD distance to the weaker 
wind star will be given by 


R2 


i + Vv~ 


:d. 


(16) 


Then, the CD will be located closer to the weaker star 
for systems with smaller r/. Furthermore, the weaker wind 
will be less diluted before the collision producing a denser 
slab (compared to systems with same r/ but different stellar 
separation). In this way, stars with large momentum fluxes 
in their winds can be important in asymmetric encounters 
as they could “force” weaker winds to radiate their energy 
rapidly, which would result in radiative shocks. In Figure 6 , 
we highlighted the stars with largest momentum fluxes in 
their winds, including their names on the plot. Encounters 
of these stars with others, of the left side of the plot, will 
produce encounters with small r/ and short distances R 2 
from the weaker star to the CD. For each of the encounters 
we registered in the previous subsection, we calculated the 
momentum flux ratio and the R 2 value. We can now modify 
the definition of the cooling parameter (Equation 3), so it 
uses the distance R 2 at which the CD is expected to form 
in the asymmetric close encounters, 


1 Vi A R 2 ,12 

Xa Sym « 2 


(17) 


where i? 2 ,i 2 = R 2 /W 12 cm. Since R 2 is always a fraction of 
d/ 2 , the density of the weaker wind will be higher at that 


position compared to the symmetric case and the slab will 
be able to radiate its energy away more rapidly. All these 
estimates are included in Table 4. From those results we 
can check that for two of our models (IDisc and Min-ecc) 
there is only one encounter where the weaker wind produces 
a cold slab that might become unstable, while in the other 
case ( 2Discs ), there are four such encounters. The difference 
is at least partially due to the fact that in the 2Discs case the 
stars are closer together and more encounters are produced 
in general. These systems deserve more study because they 
could be clump sources. It is important to remark that, in 
the asymmetric case, on one side of the CD we have a thin 
shell while on the other we have the hot shocked gas of the 
stronger wind. The latter tends to stabilise any instability 
possibly excited, so even if Xasym < 1 there might be no 
clump formation. Also notice that all these estimates were 
done under the assumption that the stars are separated well 
enough in order to accelerate their winds up to their terminal 
velocities. However, for extreme values of 77 that might not 
be the case. 


5 DISCUSSION 

5.1 Limitations and uncertainties in the model 

Although our model has been very useful to test the 
likelihood of clump formation in stellar wind shocked gas, 


MNRAS 000, 1-11 (2015) 










10 D. Calderon et al. 


there are two important assumptions we have to consider 
before comparing with other works, specially with 2- or 
3-dimensional models. 

5.1.1 The planar winds assumption 

Our model considers planar shocks, rather than more 
realistic spherical wind collisions. A planar shock geometry 
is a good approximation for the wind collision at the apex 
(i.e., at the CD and two-star axis intersection), where the 
gas effectively moves along the axis joining both stars. 
Off this axis, the winds will have a perpendicular velocity 
component, and the planar approximation breaks down. 
Moreover, the off-axis velocity will prevent the material from 
accumulating in the slab. Nevertheless, the well-modelled 
apex is arguably the most interesting location of the wind 
collision for us, as this is where we expect the most massive 
clumps to form. As shown in Figure 5, as long as we stay in 
the radiatively-efhcient regime, higher values of x generate 
more massive clumps. So far, we have defined \ using the 
wind speed and the distance between the stars, which is 
only appropriate along the two-stars axis. In a more realistic 
model, x changes along the CD, decreasing from the apex 
as we show below. 

Let us define a more general cooling parameter as 
X = xOi where 9 is the angle between the two-star axis 
and a line connecting one star and an arbitrary point P on 
the CD. The distance from the star to P will be given by 
d' t = d*/cos6*. Furthermore, the shock at P will be weaker 
than at the apex, as only the component perpendicular to 
the CD will contribute to it. This component will be given 
by V' = U cos 9. Using these quantities we can generalise 
Equation 3 obtaining 

X(0) = i^C^(cos0) 4 ' 4 = X-l(cos0) 4 ' 4 , (18) 

where x± is the cooling parameter at the apex, as previously 
defined in Equation 3. From this new expression we see 
that in general, xO ^ Xx> recovering Equation 3 when 
9 = 0 as expected. Thus, away from the apex the cooling 
parameter decreases, and according to our results, we expect 
the creation of less massive clumps. With this in mind, the 
clump masses we quote are the largest we can expect from 
a given wind collision and they would form near the apex 
where the planar-wind assumption is accurate. 

5.1.2 Impact of metallicity on radiative cooling 

Different metallicities result in different cooling timescales, 
which in turn modify the sizes and masses of the clumps. 
As there is no agreement on the metallicity measurements 
for the massive stars in the Galactic Centre (see the review 
by Genzel et al. 2010) we have set Z = 3Zq, following 
Cuadra et al. (2005) and close to the values Martins et al. 
(2007) studied. However, in order to quantify the sensitivity 
to our choice, we also produced models with metallicities 
of Zq and 5 Zq in the analytical cooling function we use 
(see Equation 10). We found changes of a factor ~ 2 in the 
clump sizes, which translate in a factor ~ 8 for the clump 
masses. As expected, a lower metallicity value increases the 
cooling timescale, and viceversa. Nevertheless, these changes 


in the metallicity do not result in switching the wind regimes 
for the system we have analysed (e.g., from radiative to 
adiabatic or in the opposite direction) as the impact on 
the cooling parameter is not very strong for the metallicity 
values we have tested. Thus, we are confident our results 
do not depend strongly on our chosen metallicity. Still, 
one should be cautious when comparing with other works 
which could have used other metallicity values, and therefore 
another cooling function. 

5.2 Binary stars 

Based on our results, clump formation seems not very likely 
to occur in the environment of the Galactic Centre, due to 
either symmetric or asymmetric stellar encounters. However, 
we have not studied a case that probably deserves more 
attention: colliding winds binary systems. 

As shown in Section 4.1, the formation of cold slabs 
typically requires stellar separations below 1 mpc. While 
those separations are not often achieved in stellar encounters 
of single stars, they are easily reached by close binaries. The 
census of young massive binary stars in the Galactic Centre 
is still incomplete, but the recent study by Pfuhl et al. 
(2014) increased to three the amount of confirmed binary 
systems:(i) IRS 16SW, a 19.5-day period Ofpe/WN9 
eclipsing contact binary (Martins et al. 2006), ( ii ) IRS 
16NE, a 224-day period Ofpe/WN9 binary and (in) E60 4 , 
a 2.3-day eclipsing contact Wolf-Rayet binary. The inferred 
separations of these binary systems are of the order of 10 ppc 
and below. That, together with the wind properties listed in 
Table 3, could make these three binaries very effective sites 
for clump formation. IRS 16SW is of particular interest, as 
it has a clockwise orbit that roughly coincided with G2’s 
at the latter’s apocentre. Notice that the wind properties 
of these stars were assigned by Cuadra et al. (2008) as 
mentioned in Section 4. However, we would not expect any 
important wind component to be substantially faster, as it 
would show up as broader lines in the spectra. Additionally, 
from the estimations by Pfuhl et al. (2014) we expect an 
overall ~ 30% spectroscopic binary fraction for the massive 
OB/WR stars in the Galactic Centre. Thus, even though our 
current knowledge of the binary population in the Galactic 
Centre is limited, close binary systems remain as a very 
promising possibility for the creation of cold clumps. A more 
detailed numerical model of this process is deferred to a 
forthcoming paper. 


6 CONCLUSIONS 

We have developed a simple and straightforward 
prescription to study clump formation through the 
NTSI mechanism in symmetric colliding wind systems. The 
input parameters are the mass loss rate, wind terminal 
velocity and stellar separation. Radiatively efficient wind 
collisions are capable of creating cold gas clumps (10 4 K) 
in a wide range of masses, where the most massive ones 
are of the order of the Earth mass for strong outflows 
(M ~ 1CP 5 Mq yr _1 ), relatively slow wind terminal 

4 The star with ID 60 in Table 3. 
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velocities (250 — 750 km s -1 ) and short stellar separations 
(0.1 — 10 mpc or 20 — 2060 AU). Nevertheless, the wide range 
of unstable wavelengths that are excited prevents us from 
predicting unequivocal clump masses, as shorter-wavelength 
perturbations grow faster and might hinder the development 
of larger-scale ones. 

We also found that the possible clump masses depend 
strongly on the timescale needed for the slab to collapse 
within the radiative wind regime. The most massive ones 
would be generated in systems where the shocked gas does 
not cool instantaneously, i.e. systems with \ approaching 
unity. Studying that regime however is not straightforward 
as both radiative and adiabatic cooling are important. 

Our results show that the formation of gas clumps 
with masses comparable to the G2 cloud is indeed 
possible in symmetric colliding winds. However, this scenario 
does not seem likely in the Galactic Centre given the 
currently known mass-losing star sample, as the required 
sub-mpc separations are very rarely achieved by them. We 
also discussed clump formation in asymmetric encounters, 
finding that the massive and slow outflow of IRS 33E could 
create clumps if confined by a powerful wind of another 
star. Similarly, the collision of winds from IRS 13E2 and 
13E4 could also generate a cold slab unstable to the NTSI. 
These stars have similar orbits and spend a significant time 
at short separations, which could explain the presence of 
many dusty clumps in their vicinity (see Fritz et al. 2010). 
However, they orbit Sgr A* in the opposite sense as G2, 
making this pair an unlikely origin for this particular cloud. 
A promising possibility is that clumps are produced in 
close binaries, of which three are currently know. Still, 
better observational data are required to constrain the wind 
properties of both components of each binary. The IRS 
16SW binary is of particular interest, as its orbit coincides 
with G2’s at apocentre. 

In conclusion, given our current analysis and the 
available stellar wind data, the formation of G2-like clouds 
in the Galactic centre appears as a possible but not very 
common event. We defer more concrete results to a future 
study using 2D and 3D numerical modelling. This is required 
to treat systems with unequal stellar winds, with x ~ 1, 
or for winds that collide before reaching their terminal 
velocity, such as compact binaries. Numerical models are 
also required to follow the growth of different unstable 
wavelengths and obtain a clump mass function. 
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